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We study the scattering resonances between two confocal hyperbolae and show 
that the spectrum is dominated by the effect of a single periodic orbit. There are 
two distinct cases depending on whether the orbit is geometric or diffractive. A gen- 
eralization of periodic orbit theory allows us to incorporate the second possibility. In 
both cases we also perform a WKB analysis. Although it is found that the semiclas- 
sical approximations work best for resonances with large energies and narrow widths, 
there is reasonable agreement even for resonances with large widths - unlike the two 
disk scatterer. We also find agreement with the next order correction to periodic 
orbit theory. 
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In recent years there has been growing interest in understanding the extent to which 
knowledge of classical mechanics can be used used to understand quantum systems [|I|. 
For sufficiently hyperbolic systems, periodic orbits provide an efficient means of calculating 
spectra semi-classically. One class of problem which has been fruitfully studied is the n- 
disk scatterer in two dimensions It has been shown that in some situations periodic 

orbit theory gives the energies and widths of the scattering resonances to great accuracy - 
typically several decimal places. 

However, this accuracy is usually only for the leading family of resonances. There exist 
other resonances with larger widths deeper in the complex momentum plane. These are badly 
approximated by periodic orbit theory because it fails to consider classical paths which, 
although not trajectories, satisfy the stationary phase condition. Recent work has 
shown that inclusion of diffractive paths - so-called creeping orbits - recovers the qualitative 
features of the exact spectrum, including the lower order families. However, the quantitative 
agreement is not as good as for the leading family of resonances which are unaffected by 
creeping. 

It is useful to consider a system which has no creeping so as to study the full spectrum 
of resonances. One such system is a pair of confocal hyperbolae. Like the two-disk scatterer 
there is only one periodic orbit, which is unstable. However, unlike the two disk scatterer, 
there are no creeping orbits. Since the system is separable we can apply WKB techniques 
as well as periodic orbit theory but this is more difficult and is not as intuitive. In addition, 
we can evaluate the resonances for the special case where one or both hyperbolae are close 
to half planes. Then we must consider edge diffiaction and this is included in the periodic 
orbit analysis. 

We want to solve the Schrodinger equation in a domain between two confocal hyperbolae, 
as shown for example in Figs. 1 and 2. This is simply the Helmholtz equation (V^ + 
/c^)^ = with the boundary conditions that on the two hyperbolae vanishes and that 
in the region between them \1' approaches /(0) exp{ikr) / y/r for large r (where r and are 
polar coordinates.) We will work in hyperbolic-elliptic coordinates, fi and 6, defined by 
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X = a cosh fi COS 6 and y = asinh/isin^ 0. Curves of constant /i are ellipses and curves of 
constant 6 are hyperbolae. In both cases the foci are at y = and x = ±a. We take the 
coordinates to be in the ranges < < vr and — oo < /i < oo and label the right and left 
hyperbolae respectively by 6i and 62 so that 62 > 61. We define length units such that the 
closest distance between the hyperbolae is 1. It follows that a = l/(cos6'i — cos 6*2). This 
can be generalized to any spacing L by substituting kL for k in what follows. 

Using the separation \E'(r) = M{fi)Q{6) and the expression = (^ + ^) /a^(cosh^ /i— 
cos^ 9), we obtain the Mathieu equations 



d^M 



+ a\b' - k' cos' e)Q = (la) 



+ a^(-6^ + rcosh>)M = 0, (lb) 



where a^6^ is a separation constant. Although we can work with these equations, it is 
convenient to express equation (|bD as a one dimensional potential problem. This is done 
by a change of coordinates ( = asinh/x and M^jj) = Z{Q/{('' + a^Y^^ yielding 

d^Z 



with an effective potential 



+ {k' - 2V{0)Z = 0, (2) 



The asymptotic boundary condition is Z{Q ~ ex^i^^O- There is a y — —y symmetry so 
there are two symmetry classes; Z'{Q) = for even states and Z{0) = for odd states. The 
effective potential is repulsive so there are no bound states, only resonances. 

Resonances occur for negative imaginary k which means that asymptotically the wave 
function is increasing exponentially [0 as 0{exp{kiQ) where k = k^—iki. Such wavefunctions 
are called Siegert or Gamow states and their asymptotic nature is a reflection of the time- 
dependent decay. To work with such wave functions numerically, one must complexify the 
coordinate ( = \(\ exp(ia) [|]. Asymptotically the wavefunctions are decaying exponentially 
function of |C| if a > — arg(A;). 



We numerically solve the simultaneous equations for k and b by shooting 0. The exact 
spectra for two choices of parameters are shown in Figs. 1 and 2 as open diamonds. Selected 
data are also listed in Tables I and II. In Fig. 1, 6i = 0.4 and 62 = 0.9 which is a smooth 
system since the radii of curvature of the hyperbolae are comparable to the inter-hyperbola 
spacing of 1. In Fig. 2, 61 = and 62 = 2.0 which is sharp since the right hyperbola is a half 
plane and its radius of curvature is zero. As will be discussed, the theoretical understanding 
of the two cases is quite different. For each value of k there is one corresponding value of b. 
The spectra of k and b are similar in both cases; their real components are very close and 
their imaginary components differ approximately by a constant. 

Periodic orbit theory works by approximating the trace of the Green function, the poles 
of which are bound states or scattering resonances [||]. This has been worked out for the 
two disk problem in terms of the single periodic orbit in the system. The result, adapted 
to this system and including both symmetry classes, is 

Kj = {n + 1)71 - i ' ^ log A + ■ ■ ■ (4) 

where n, j = 0, 1, • ■ ■ and A = tan^(^^2/2)/ tan^(0i/2) is the stability for a complete traversal 
|n|. The y —y parity is (—1)-'. 

Equations ( |Ta|) and (^) can also be solved using WKB theory. The details will be pre- 



sented elsewhere |10]. We begin with the expansion 



k = K-iki + 0{k-^) (5) 
b = br-ibi + 0{b-^). 

b is the eigenvalue of the equation (£+6^)\I' = where C = (cosh^ /i^+cos^ / (cosh^ /x— 
cos^ 9) as can be seen from equation (1). In the neighbourhood of the periodic orbit and to 
leading order, the equation for b is the same as the equation for k so that b,. = kr- 



The WKB approximation for ([T^) is [|ri|] 



(n + l)7r = a / dO^b^ - k^ cos^ 9 (6) 



This equation can be expressed in terms of incomplete elliptic integrals [|T2| . We expand the 
integrand of equation (|^) in powers of (6^ — k^) / (/c^ sin^ 9) and keep the first two terms. This 
gives kr = {n + l)7r and bi = fki where / = (1 — 2/alogA). The imaginary components of 
k and b are related by a factor / which depends only on the geometry and is the same for 
all resonances. 



The WKB approximation for the resonances of equation (y) is |Tl 



(2j + l)7r = 2a / dzJk'^ - 2V{z) (7) 

J Z- 

where V{z) is shown in equation (|^). This formula comes from identifying the resonances as 
the poles of the transmission coefficient. Zj- are complex turning points which are solutions 
of V{z±) = k"^ /2. We arrive at a solution in which k is as in equation and 

bnj = (n + l)7r - tf'^^^ log A + ■ ■ ■ (8) 

These results are consistent with periodic orbit theory. However the calculation of b is 
new. The imaginary component of b is less negative than that of k while the real components 
are equal to leading order. The semiclassical results are shown in Fig. 1 as crosses and are 
also listed in Table I. The worsening of the agreement as we go down in the k plane is 
expected [|13| and is consistent with the expansion (|). 

It is interesting to study the higher order terms in the expansion of the resonances. The 
first two corrections to the leading family of resonances yield [jl4|,|15 



[n + IjTT \ 4 (n + lj^vr^y 



where Cr and q are constants. For the symmetric case 62 = ti — 6i, v^e have [|1^] = cot^ 61. 
For the example 61 = 1.0 we fit the 200'th resonance to equation to extract c,. = 0.4123 in 
agreement with the expected result. We also find Cj = 0.593 but there is as yet no analytical 
result. The two constants have been worked out for the two disk problem |T^ where it is 
found that q is positive. This means that the resonance lifetimes are longer than the first 
order approximation, as is also found here. Since periodic orbit theory is an asymptotic 



expansion in h we expect that the higher order corrections will eventually blow up for 
any resonance n. The results for the scaling with j are ambiguous so we leave this for a 
later publication. 

Note that some of the resonances are not in the physical region kr > but instead have 
kr < 0. (Because the dimension is even, the pole structure is not symmetric with respect 
to reflections through the imaginary k axis ^^■) For example, the state with indices n = 
and j = 4 does not exist as a physical resonance although there is a semiclassical prediction. 
More resonances lie in the unphysical region for large j and for small 6i. 

For the configuration of Fig. 2 we need a different analysis since equation predicts 
infinite widths. The previous results relied on an expansion in powers of 1/ sin^ 6' which 
breaks down if ^ 1 or tt — ^ 1- In that case we can approximate the hyperbola by a 
wedge so that a qualitatively new phenomenon is responsible for the resonances. There is a 



class of orbits which are not classical trajectories but rather are diffractive [|T8| as shown in 
Fig. 2. The ray starting at P illuminates the vertex which then acts as a point source for an 
outgoing circular wave which illuminates P'. The Green function for this process ]TB|,|5| is 



G(P',P,fc) = -^«4±^G(KP.*) (10) 

y/Zirkr 

where G{V, P, k) is the free space Green function which connects the point P to the vertex 
V and can be approximated by the standard Van Vleck formula. The geometric factor d 
depends on the incoming and outgoing angles and on the wedge angle. 

As with geometric orbits taking the trace of G selects the closed orbits. These are all 
repetitions of the primitive orbit which diffracts from the wedge and bounces back off the 
wall, with a sign change. When the left wall is straight the primitive orbit contributes to 
the trace an amount —igQ{k)/A, where 



go{k) = d exp(z {2k + tt/A)) /V4^. (11) 



For normal incidence d = Y^^j^^^^yfi/l with 7 = 2(1 — ^^i/vr). Note that d = 1 for 
6*1 = 0. We must also include the focusing properties of the left wall. The theory of reference 



18| relies on a picture of cones such that amphtudes decrease with the square root of the 
cone widths. For a straight wall, a cone leaving the vertex with an opening angle (f) returns 
with a width 20. For a curved wall, the cone returns with a width 20(1 + where R = 
— sin 6*2 tan 6*2/(1 — COS62) is the radius of curvature of the left wall. Therefore, the amplitude 



is divided by a factor of Jl + 1/R and the denominator of equation (O) is replaced by 



47rA;(l + 1/R)- This factor can also be derived from the integration perpendicular to the 
diffractive orbit when the trace is evaluated but the geometric argument above is clearer. 
Multiple diffractive terms in the trace are just powers of the primitive orbit ^ so 

• 00 

TrG(fc) = -7E^?S(^) (12) 

^ r=l 

i/A 



exp{-i{2k + 7r/4) + log ^ V47rA;(l + 1/R)} - 1 
The poles are given to leading order by 

7, z, 47r2(n + 7/8)(l + l/i?) 
kn = {n + -)n - - log ^ '-^ ^ (13) 

The factor of 7/8 comes from the n/A phase shift in equation (^) and is an important 
qualitative distinction from the geometric theory. A simple extension to the case where 
both hyperbolae are sharp yields 

3, z, 27r2(n + 3/4) 
fc. = (n + -).-- log (14) 

A WKB analysis is possible for the sharp geometry. For simplicity, we just discuss the 
result for 9i = 0. Then the integral in equation (H) can be expressed as a complete elliptic 
integral of the second kind ||12| and expanded in powers and logarithms of (6^ — A;^)/6^. 
Again we use the expansion (|^) and argue that = kj.. The WKB analysis of equation 
(0) is the same as before. The result of these calculations is that k is given as in equation 
(p^) and bn = k„ + i/2a. As with the geometric case, the imaginary component of b is less 
negative than the imaginary component of k. 
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The results of the semiclassical analysis for the diffr active case are shown in Fig. 2 and 
in Table II. We find good agreement with the exact results for both k and b. As in the 
geometric case, the exact lifetimes are longer than the lowest order semiclassical prediction. 
The theory only predicts one family of resonances. This is consistent with the numerics 
where we see the that widths of the next family increase logarithmically as ^ 0. 

There are two unexplained features of the spectrum for 9i small but nonzero. One is 
the crossover to the geometric result near ~ 1/{t^9\) (by the criterion that 47r|/c| ~ A). 
The other is the positions of the lower families. These might be explained by including both 
geometric and diffractive orbits in the Green function. These calculations should be possible 
for the WKB analysis as well. Another interesting issue is the higher order corrections which 
are calculable using either periodic orbit theory or WKB. We can also use periodic orbit 
theory to approximate the resonance spectrum in a system comprised of wedges of finite 
angles. These issues will be addressed in future publications. 

I would hke to thank Gabor Vattay for suggesting this project and for showing me the 
results of his unpublished calculations. I benefitted tremendously from discussions with him 
and with Stephen Creagh and Andreas Wirzba. I would also like to thank Charles Clark for 
telling me about the complex coordinate rotation method. This work was supported under 
the EU Human Capital and Mobility Programme. 
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FIGURES 

FIG. 1. The top box is the configuration space for 6i = 0.4 and 62 = 0.9. The sole periodic orbit 
is shown as a dotted Une. The middle and bottom boxes show the spectra of k and b respectively, 
plotted in the complex plane. The exact quantum resonances are represented as open diamonds 
and the semiclassical predictions are represented as crosses. 

FIG. 2. The same as Fig. 1 for 61 = and 62 = 2.0. The dotted line in the top box is the 
diffractive periodic orbit. The crooked dashed line is a typical diffractive path connecting the 
points P and P' . Its Green function is given by equation (|lO|). 
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TABLES 

TABLE I. The exact and semiclassical predictions of k for a selection of resonances and for the 
geometry of Fig. 1. 

n j k k^'^ 

3.0587 -zO.3930 3.1416 -zO.4342 

2 9.3881 -iO.4258 9.4248 -iO.4342 

10 34.5467 -iO.4334 34.5575 -iO.4342 

20 65.9677 -i0.4340 65.9734 -iO. 4342 

6 21.9714 -iO.4326 21.9911 -iO.4342 

6 6 21.2799 -i5.5685 21.9911 -z5.6443 

6 12 19.4438 -zlO.4733 21.9911 -zlO.8544 



TABLE IL As in Table I for the geometry of Fig. 2. 

n j k k^'^ 

2.6537 -il.0060 2.7489 -il.0201 

2 8.9913 -il.3160 9.0321 -zl.3175 

10 34.1514 -il.6501 34.1648 -^.6501 

20 65.5731 -il.8131 65.5807 -il.8131 
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